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ABSTRACT 

We report a spatially resolved, X-ray spectral analysis of NGC 4258 using 
archival Chandra and XMM-Newton observations. The XMM-Newton spectra 
of the nuclear region are well described by two power-law components, a soft 
(0.57 keV) thermal component, and an Fe Ka line with EW = 40 ± 33 eV. The 
properties of the second, weaker power-law component are similar to those of an 
off-nuclear source 2.5" SW of the nucleus. The spectrum of the extended emis- 
sion of the entire galaxy is well described by two thermal components (MEKAL) 
models with temperatures ~ 0.60 and 0.22 keV. The Chandra and XMM-Newton 
spectra along the anomalous arms show that the absorbing column density to the 
SE anomalous arm is consistent with absorption by gas in our Galaxy, while the 
absorbing column to the NW anomalous arm is higher, indicating that the NW 
arm is partially on the far side of the galactic disk. The combined Chandra data 
clearly detect the X-ray emission from the hot spots at the end of the approxi- 
mately N-S radio jets. By assuming the hot spots represent shocked thermal gas 
at the ends of the jets, we estimate shock powers of ~ 3 x lO 39 /" 1 / 2 erg s _1 (/ 
is the filling factor), similar to the radiative power in the inner anomalous arms, 
consistent with the notion that the jets could be responsible for heating the gas 
in the anomalous arms. 

Subject headings: galaxies: individual: NGC 4258 — galaxies: Seyfert — galaxies:spiral— 
X-rays: galaxies 

1. Introduction 

The interaction between nuclear activity and host galaxy is a very interesting problem 
the study of active galactic nuclei (AGNs). Such interactions could be rather strong as 
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in the case of NGC 4258 (M 106), a nearby (7.2 Mpc, Herrnstein et al. 1999) prominent 
barred Seyfert 1.9 galaxy. NGC 4258 is one of the best laboratories to study accretion onto 
blackholes because of its well determined blackhole mass and distance. The spatial-velocity 
distribution of water mega-maser sources in NGC 4258 on scales of 0.14-0.28 pc indicates a 
thin Keplerian disk rotating around a putative blackhole with mass 3.9 x 10 7 M Q (Miyoshi 
et al. 1995; Herrnstein et al. 1999). The accretion disk has an almost edge on orientation 
with the rotation axis pointing almost north (PA = -8°) and tipped 82° to the line-of-sight . 
A pc scale jet closely aligns in projection on the sky with the rotation axis (Herrnstein et al. 
1997). 

The galaxy is well known for its pair of "anomalous arms" , located between the normal 
spiral arms of the galaxy, and visible in Ha, soft X-ray, and radio continuum. These arms 
emerge from the nuclear region as linear features, then bend at 2-5 kpc radius in a trailing 
sense with respect to the galactic rotation. The smooth appearance and the lack of blue 
light from young stars make them distinctive from the normal spiral arms. Many authors 
suggest that the anomalous arms are ejected from the nucleus (e. g. van der Kruit, Oort, & 
Mathewson 1972; van Albada 1977; Icke 1979), in the forms of either "cannon balls" or jets. 
However, there is strong evidence that at least the inner part of the anomalous arms lie in the 
disk of NGC 4258 (van der Kruit 1974; van Albada & van der Hulst 1982; Hummel, Krause, 
& Lesch 1989; Plante et al. 1991; Cecil, Wilson, & Tully 1992). It is hard to understand 
how the ejecta could be confined to the thin galactic disk while propagating to the observed 
distances of > 2 kpc. The energy requirement to push aside the dense interstellar gas in the 
disk is also unrealistically high. A more serious problem with the jet models is that a pair 
of currently active jets, which are extensions of the parsec scale jet seen in VLBA images 
and which extend to projected distances of 870 pc south and 1.7 kpc north of the nucleus, 
is clearly visible in the VLA image (Cecil et al. 2000). The galaxy disk has an inclination 
of 64°, while the radio jet lies very close to the polar axis of the maser disk, which itself is 
inclined by 60° to the polar axis of the galaxy stellar disk. This misalignment between the 
currently active jets and the anomalous arms makes the jet models for the anomalous arms 
unlikely. X-ray emission from the anomalous arms was first discovered with the Einstein 
HRI (Cecil et al. 1992). Subsequent observations from ROSAT (Pietsch et al. 1994; Cecil, 
Morse, & Veilleux 1995; Vogler & Pietsch 1999) have found that the X-ray emission is 
closely associated with the anomalous arms, and is thermal in nature. The temperature of 
the hot gas was found to be kT ~ 0.3 — 0.6 keV, with a possible halo component (Vogler & 
Pietsch 1999). Combining the results from VLA and VLBA, and with the help of a 14 ks 
Chandra observation, Wilson, Yang, & Cecil (2001, hereafter WYC) showed that the vertical 
projection of the jets onto the galactic disk aligns closely with the inner anomalous arms, 
suggesting the anomalous arms may be shock heated gas powered by hot halo gas driven into 
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the galactic disk by the out-of-plane jets. To better understand the nature of the anomalous 
arms, their relation with the jets and the AGN, we need detailed knowledge of the X-ray 
emitting gas in the anomalous arms, the jets and the nucleus. Such knowledge will allow 
constraints on the density, pressure, and energetics in the jets and the anomalous arms. 

The nucleus of NGC 4258 is heavily obscured in the ROSAT band, and the first hard 
X-ray detection was reported using ASCA observations, with the absorption column density 
Nh — 1.5 x 10 23 cm -2 . A narrow 6.4 keV Fe Ka line was detected with an equivalent 
width (EW) of 250 ± 100 eV (Makishima et al. 1994). Fruscione et al. (2005) noticed that 
the column density remained mostly constant (~ 10 23 cm' 2 ) over 9 years, except on two 
occasions, which they attribute to absorption by a warped accretion disk. On the other 
hand, flux variability in the 2-10 keV X-ray band has been seen on timescales from 19 days 
to years, with little change in the power-law index (Reynolds, Nowak, & Maloney 2000; 
Terashima et al. 2002; Fruscione et al. 2005). The Fe Ka line was found to be weaker in 
recently reported ASCA (Reynolds et al. 2000), BeppoSax (Fiore et al. 2001) and RXTE 
(Mushotzky, private communication) observations, with EW~ 100 eV. Recent observations 
from XMM tend to find much lower EW, and in many cases, only upper limits (Pietsch & 
Read 2002; Fruscione et al. 2005). Moreover, weak Fe Ka absorption features were reported 
using a Chandra observation (Young & Wilson 2004), but not detected in XMM-Newton 
observations. The variability of EW has been confirmed in a recent reanalysis of the archival 
ASCA, RXTE and XMM data, using the latest calibrations (Hanson et al. 2006, submitted; 
hereafter Hanson06). 

In this paper, we present our analysis of the archival Chandra and XMM-Newton ob- 
servations of NGC 4258, with a focus on the anomalous arms and jets. With the excellent 
spatial resolution of Chandra and the sensitivity of XMM-Newton, we are able to probe X-ray 
emitting gas in the anomalous arms and jet hot spots in detail. In § we first present the 
data reduction. The results are presented in § El We discuss the implications of the results 
in § HI and give concluding remarks in § [5l 

2. Data reduction 

2.1. Chandra observations 

NGC 4258 was observed with the Advanced CCD Image Spectrometer (ACIS) on Chan- 
dra on April 17, 2000 and May 28, 2001 (Tabled]). The nucleus was located at the aim point 
on the S3 chip. The spatial resolution in the 0.4-2 keV band is < 1" for all the off-axis angles 
covered by the detected diffuse emission from NGC 4258. The default frame time of 3.2 s 
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was used. The data were reduced using CIAO 3.3 with calibration version CALDB 3.2. We 
followed the standard procedure in reducing the data. By comparing the positions of the 
nucleus and point sources which have radio counterparts in the 6 cm high resolution Very 
Large Array (VLA) image used in Cecil et al. (2000), we found the astrometric error of the 
Chandra image is < 0.5". 

In Fig. [TJ we show the adaptively smoothed "true color" image of NGC 4258 in the 
0.4 - 2 keV band. The anomalous arms appears smooth but not featureless. The northern 
hot spot (marked "N") of the radio jet is clearly visible in the X-ray image. An arc-like 
structure, suggestive of a shock wave, is seen ~ 6" to the west of the northern hot spot 
(N) in the northern anomalous arm, indicating a possible relation with the hot spot. The 
arc is ~ 1' from the nucleus. A similar "arc" is also visible, but less prominent, in the 
southern anomalous arm, at about the same distance from the nucleus as the northern arc. 
The X-ray color is not uniform along the anomalous arms and the surrounding regions. It is 
obvious that the SE arm is in general redder (softer) than the NW arm. The diffuse emission 
surrounding the anomalous arms is generally soft, with the emission from the NE quadrant 
redder (softer) than that in the SW quadrant. 

We extract spectra in regions along the anomalous arms as well as the hot spots which 
are detected in X-rays (Fig. [2]). The regions are chosen so that there is little variation of color 
within each region. In the presence of large scale X-ray emission, extra care is needed in 
determining the level of background. We compare the following methods. First, we visually 
inspect the image and extract the background spectrum from two regions with no point 
sources and which show little diffuse emission. The source spectrum from a region on the 
anomalous arms with low surface brightness is also extracted. The resulting two background 
subtracted spectra using the two background files are compared, and no significant difference 
was found. Second, we extract background spectra at the same physical location as the 
source from a stack of background files compiled by Maxim Markevitch0 , which has been 
reprojected to the same sky coordinates as our observations. The resulting background 
subtracted spectra is compatible with those using the first method. 



2.2. XMM-Newton Observations 

We use the data from five XMM-Newton observations of NGC 4258 available in the 
archive (Tabled]). The data is reprocessed using the latest calibration files and SAS 6.5.0. 
The resulting event lists from the two MOS detectors and the PN detector are filtered for 
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background flares, using visual inspection of the light curves of the whole field. The resulting 
good exposure times can be found in Table [TJ Comparing with the filtered exposure time 
of Fruscione et al. (2005), which uses a signal-to- noise criterion, we found good agreements. 
The total good time from the observations are 132.4 ks for the two MOS detectors combined 
and 45.3 ks for the PN detector. 

We extract spectra along the anomalous arms and from the regions shown in Fig. [2j 
which are determined using the Chandra data (left image). The background spectra are 
taken from regions of low emission on the same chip close to the "source" regions. The 
source region sizes are large enough that the number of events from neighboring regions due 
to the point spread function (PSF) wings is small, except in regions close to the nucleus 
(particularly regions 10 and 11), where scattered light from the bright nucleus affects the 
spectra of the diffuse emission. However, this "mixing" does not seem to alter the qualitative 
conclusions for these regions (see § 13.21) . The effect can be modeled using xmmpsf in XSPEC. 
Unfortunately, the PI spectra produced by the current version of SAS are not compatible 
with the xmmpsf model. XSELECT included in the ftools package can produce spectra 
that are compatible with xmmpsf, but only spectra from a single instrument from a single 
observation can be treated properly with xmmpsf. Thus we rely more on comparing the 
XMM-Newton spectra with those from Chandra. 

The spectra from different observations which use the same detector (MOS or PN) are 
co-added using the addspec task in the ftools package. The routine combines the spectra from 
the source, background and produces a new response matrix. The spectra from both MOS 
detectors are also added. The addition of the spectra of the extended emission is justified 
because the emission is not variable. 

While the cross calibration between PN and MOS has improved significantly since the 
release of SAS 6.5.0, the calibration is still not very good below 0.8 keV and above 7 keV0. 
This systematic uncertainty can cause spectral fitting to produce larger \ 2 values if the 
spectra from both PN and MOS are fitted simultaneously. Experimenting with PN and 
MOS spectra, we found that the MOS spectra tend to produce more artifacts that mimic 
emission lines in the 0.2-2 keV band. In the case of thermal gas with temperature ~ 0.5 keV, 
fitting MOS spectra produces higher abundance values. In this study, we present the spectral 
parameters using only PN as well as using both PN and MOS detectors. 



2 see Kirsh (2006), "EPIC status of calibration and data analysis" (CAL-TN-0018). 
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3. Results 

3.1. The nuclear region 

The spectrum of the nuclear region has previously been reported using XMM-Newton 
(Pietsch & Read 2002; Fruscione et al. 2005, Hanson06) and Chandra (Young & Wilson 
2004; Fruscione et al. 2005). These observations, as well as the earlier ASCA observations 
(Reynolds et al. 2000), have produced two puzzles: 1) A weak second power-law component 
is needed in addition to the heavily absorbed power-law. The nature of this component 
is unknown. It could come from scattered light from the nucleus, or from unresolved point 
sources; 2) Weak absorption feature were identified in Chandra observations but not detected 
in XMM-Newton observations (Fruscione et al. 2005 reject these detections, but the issue 
seems to lie in the method of statistics). 

As noticed in WYC and Young & Wilson (2004), a bright off-nuclear source (J121857.3+471812) 
2.5" to the SW of the nucleus has been found in the Chandra observations. Because the source 
is significantly weaker than the nucleus, it is commonly assumed that it contributes little 
to the spectra obtained with instruments which cannot separate the point source from the 
nucleus, and it is thus ignored. In this section, we re-examine the spectra of the nuclear 
region and the off-nuclear source, and discuss the effects of the off-nuclear source on the 
XMM-Newton spectrum in which the off-nuclear source is not resolved from the nucleus. We 
then compare the Chandra spectra of the off-nuclear source and that of the nuclear region 
from XMM-Newton. In contrast to previous works, we use stacked spectra to improve the 
signal-to- noise ratio. 

3.1.1. The XMM-Newton Spectra of the nuclear region 

The spectrum in each observation is extracted from a circular region centered at the 
nucleus with a radius of 20". To obtain meaningful results from a stacked spectrum, either 
the source or the response of the instrument has to be invariable. While little variability is 
seen below 2 keV, the spectra of the nucleus show moderate variability above this energy. We 
check the variation of the EPIC responses using simulations. The response matrices produced 
for each of the observations are convolved with the best-fit model from the longest XMM- 
Newton observation. Above 2 keV, there is no difference between the convolved spectra, 
indicating little change in the response in the hard band. Below 1 keV, on the other hand, 
there is a small change in the response. Fortunately, the spectra in this band show no 
variability. This validates the stacking method. 
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We first study the stacked PN spectrum in the 0.2-10 keV band. The spectrum is 
grouped to a minimum of 25 counts per bin to allow use of \ 2 fitting. We employ Model 
J from Hanson06, which is composed of two power-law components, a collisionally ionized 
thermal component based on the MEKAL model (Mewe, Gronenschild, & van den Oord 1985; 
Mewe, Lemen, & van den Oord 1986; Liedahl, Osterheld, & Goldstein 1995), and an Fe Ka 
emission line. The first power-law component comes from the AGN and is more obscured. 
The second power-law is much weaker than the first one, and is assumed to be unobscured 
by the gas shrouding the AGN. The best-fit parameters and the 90% confidence ranges for a 
single interesting parameter (in small numbers next to the parameters) are listed in Table [2] 
(Fit 1). The aperture corrected model 2-8 keV flux is (7.1 ± 0.4) x 10~ 12 ergcm _2 s _1 (la 
error). The corresponding unabsorbed luminosity is 4.4 x 10 40 ergs -1 . 

To better constrain the Fe Ka emission, we use both the PN and MOS spectra above 
2.2 keV (Fig. [3]). Choosing only the hard band emission allows us to use a single absorbed 
power-law as the baseline model of the continuum. We obtain Nh = (8.2 ±0.3) x 10 22 cm~ 2 , 
T = 1.57 ± 0.05 §. A narrow iron Ka emission line (Eu ne = 6.39 ± 0.04 keV, cr < 0.1 keV) is 
clearly seen in both the PN and MOS spectra. By adding a narrow emission line (3 degrees 
of freedom, dof.) \ 2 improves from 992.5 to 969.2, corresponding to a confidence level of 
99.996% according to an F-test. The equivalent width of the line is 40 ±33 eV. The line flux 

IS 3.^2 2 ^ 

10 6 photons cm 2 s 1 . These values agree very well with those from Hanson06, 
in which the PN data from the five observations have been fitted simultaneously, but with 
the absorption left free for each spectrum. The properties of the Fe Ka emission line also 
agrees with the upper limit found in Fruscione et al. (2005). 



3.1.2. Off-nuclear source 

The off-nuclear source J121857.3+471812 was first reported in WYC based on the Chan- 
dra April 17, 2000 observation. The Chandra spectrum from the May 28, 2001 observation 
has also been reported by Young & Wilson (2004). We extract the Chandra spectra of the 
source from a circular region with radius of ~ 1.5" (which corresponds to an enclosed energy 
of 95% at 1.4 keV ). The spectra of the background have been extracted from a region at 
the same distance from the nucleus but on the opposite side of the off-nuclear source. This 
allows the PSF scattered nuclear photons in the source extraction region to be properly 
taken into account. Spectra from both observations can be described very well with a single 



3 The errors represent the 90% confidence ranges for a single interesting parameter. This convention is 
used throughout this paper unless noted otherwise. 
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absorbed power-law, and the resulting parameters agree well with the previously reported 
values. There is no significant variability between the two epochs. This allows us to fit the 
combined spectra to improve the signal-to-noise ratio. The resulting best-fit parameters are 
N H = 3.4|° x 10 21 cm" 2 , T = 1.74 2 ;° 2 , / . 5 _ 8 keV = 1.2j;| x 10~ 13 erg cm" 2 s" 1 . If the source 
is in NGC 4258, its 0.5 - 8 keV luminosity is ~ 7 x 10 38 ergs" 1 . 

We compare the flux from the off-nuclear source with that of the second power-law 
component in the XMM-Newton spectra. We replace the second power-law in FIT 1 (Ta- 
ble [2]) with an absorbed power-law with the same index and absorption column density 
as those of the off-nuclear source, found from the Chandra observations. The resulting 
best-fit parameters are shown in Table [2] as FIT 2. The x 2 increases slightly in this new 
fit but is still acceptable. The new fit is still better than that without introducing the 
second power-law component. The 0.5-8 keV flux of the second power-law component is 
~ 1.4^o x 10~ 13 ergcm _2 s _1 , which agrees very well with the flux of the off-nuclear source 
from Chandra. This suggests that the off-nuclear source is likely to be responsible for the sec- 
ond power-law component in the XMM-Newton spectrum and leaves less room for scattered 
light from the nucleus. 

We also search for emission lines in the combined spectra of the off-nuclear source. Given 
the small number of photons in these observations, we regroup the spectrum to > 3 counts 
per bin. Interestingly, an emission line like feature at ~ 6.4 keV is seen in this rather noisy 
spectrum (Fig. Unfortunately, only the May 28, 2001 observation has enough counts 
above 5 keV to show this feature. By adding an unresolved Fe Ka emission line to the 
continuum model, and using the C-statistic (Cash 1979), we found the energy of the line 
feature to be Eu ne = 6.36g 26 keV (la confidence range) and an upper limit to the line flux 
of 1.6 x 10~ 6 photons cm" 2 s _1 . To check if the line could be caused by X-ray photons 
from the nucleus scattered by the HRMA, we examine the stacked Chandra image extracted 
between 6.2 and 6.6 keV. A small number of events can be seen around the nucleus. Within 
the spectrum extraction region of the off-nuclear source, 3 events are found, two of which 
are in the same pixel. In an annular region with rj„ = 1", r out = 4" and which is centered on 
the nucleus and overlaps the off-nuclear source extraction circle, we found 7 counts in the 
6.2 to 6.6 keV band, all of which are in the SW quadrant with respect to the nucleus (i.e. 
in the direction of the off- nuclear source). Comparing the counts in the annulus and the 
source extraction region, we found the significance level of detection of the off-nuclear source 
is ~ 84% between 6.2 and 6.6 keV (we have used the best-fit upper-limits in Gehrels (1986) 
to calculate the significance level). We therefore cannot rule out the possibility that the line 
represents HRMA scattered nuclear light. However, given that no strong Fe Ka emission 
line is detected in the nuclear spectrum in the May 28, 2001 Chandra observation, and that 
the photons detected around the line energy seem to "tilt" toward the off-nuclear source, it 
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is very possible that the emission line does come from the off-nuclear source. It is noticeable 
that the upper limit to its line flux from the Chandra spectrum is only a factor of 2 lower 
than the Fe Ka emission line flux of the entire nuclear region from XMM-Newton. Should 
the line be confirmed to come from the off-nuclear source, it would certainly complicate the 
interpretation of the Fe Ka line results from instruments that cannot resolve the off-nuclear 
source from the nucleus. An emission line from the off-nuclear source might also "fill-in" one 
of the putative absorption lines from the nucleus, such as those reported in Young & Wilson 
(2004), making the absorption line harder to detect with XMM-Newton. Unfortunately, Fe 
Ka absorption from the nucleus was only reported in only the April 17, 2000 observation 
(Young & Wilson 2004). This prevents a direct check on how the tentative emission line 
from the off-nuclear source would affect the spectrum of the nucleus if the off-nuclear source 
were not resolved from the nucleus, as is the case in the XMM-Newton observations. 

3.2. The Hot gas in NGC 4258: Anomalous arms and other diffuse 

components 

As shown by ROSAT observations (Cecil et al. 1995; Vogler & Pietsch 1999), the 
extended emission in NGC 4258 is predominantly thermal. Vogler & Pietsch (1999) found 
that the emission can be decomposed into two thermal components, with temperatures of 
~ 0.5 and 0.2 keV. The absorption corrected total luminosity L 01 _ 2 a kev ~ 2 x 10 40 ergs -1 . 
The low temperature component is suggested to originate from a halo component, similar 
to those found in other spiral galaxies. 

We first look at the global spectrum of the galaxy. We extract the XMM-Newton spectra 
from an elliptical region centered on the nucleus, with semi-major and minor axes of 2.6' 
and 1.5', respectively. The major axis of the ellipse is in the direction of the anomalous 
arms and is close to the major axis of the galaxy disk. The extraction region covers most 
of the soft emission while avoiding any bright point sources. Because the low energy tail of 
the absorbed power-law emission from the nucleus blends with the thermal emission around 
2 keV, we model the spectra over 0.2-10 keV. We first test a model similar to the one we 
used to fit the nuclear spectra, which consists of two power-laws, a MEKAL, and an Fe Ka 
emission line. The model produces a very good fit to the spectrum above 2 keV, and agrees 
very well with the best-fit to the nucleus. Below 2 keV, a large residual from this model is 
clearly seen, indicating the emission cannot be well described with a single thermal emission 
component. We add a second MEKAL component. Unfortunately, the abundance of this 
component is poorly constrained by the data, and the best-fit abundance is unrealistically 
high. We thus choose a fixed abundance equal to 0.15 of the solar value, a typical value 
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found in the outer regions of the extended emission (see below) . The addition of the second 
MEKAL component significantly improve the fit. The temperatures of the best-fit MEKAL 
models are 0.595 ± 0.005 and 0.22 ± 0.01 keV. This confirms the results of Vogler & Pietsch 
(1999). The 0.2-2.4 keV fluxes are found to be (2.38 ±0.05) x 10" 12 erg cm" 2 s" 1 for the 0.59 
keV component , and (5.10±0.06) x 10 -13 for the 0.22 keV component. The total absorption 
corrected luminosity is Lo. 1-2.4 keV ~ 1-7 x 10 40 ergs -1 . 

The spectra in the regions shown in Fig. [2] are investigated using both Chandra and 
XMM- Newton data. We first model the extended emission with a single thermal component, 
which produces a very good fit to all the Chandra spectra between 0.4 and 2 keV. In all 
regions except 6, 10, and 11, the XMM-Newton spectra in the 0.2-2 keV band are used. In 
these three regions, the hard band emission from the AGN contributes significantly to the 
0.2-2 keV spectra. We thus model the spectra in regions 6, 10, and 11 between 0.2 and 
10.0 keV using a MEKAL plus an absorbed power-law model. In Table |3j we list the best-fit 
parameters using XMM-Newton spectra. While producing consistent results, including the 
MOS spectra seems to systematically decrease the best-fit absorption column density and 
increase the temperature. The reduced \ 2 P er dof is also higher if the MOS spectra are 
included. This is still true even if we allow the normalization of the models to be different 
for the two instruments. 

The best-fit absorption column densities and temperatures of the thermal gas along the 
anomalous arms are shown in Fig. [5J It is obvious that Nh increases along the anomalous 
arms from SE to NW. The best-fit column densities in the SE arm are consistent with that 
of the Milky Way Galaxy in the direction of NGC 4258, indicating little absorption intrinsic 
to NGC 4258 in the SE arm. The column densities in the NW arm, on the other hand, are 
a factor of 3-4 times the Galactic value. The temperature, on the other hand, peaks at the 
nucleus (Region 6) and decreases slowly as a function of distance from the nucleus in both 
NW and SE directions. There is no systematic difference in the temperature distribution 
between SE and NW arms. The increase of column density is naturally interpreted if the 
NW arm is slightly behind the galactic disk of NGC 4258, as argued by WYC. The SE arm, 
on the other hand, is either in the galactic plane, or in front of it. 

The emissions in regions 10 and 11, which are located to the NE and SW of the nucleus, 
are contaminated by scattered nuclear light from the broad PSF wing of XMM-Newton. 
This is much less of a problem for Chandra. We first examine the XMM-Newton spectra 
of Regions 10 and 11 without using the xmmpsf model. There is little difference in the 
best-fit temperature between the two regions (Table [3]). On the other hand, the absorption 
column density in Region 11 is about twice the value of Region 10. Joint analysis of the 
spectra suggest the difference in the column density is at > 4a level. Given the inclusion of 
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scattered nuclear light in both spectra, the difference in column density is likely to be higher. 
Modeling 0.4-2 keV Chandra spectra of the two regions yields Nh = 1-7qq x 10 20 cm -2 , 
kT = 0.49°;| keV for Region 10 and N H = 7.8* 4 5 x 10 20 cm" 2 , kT = 0.57°;^ keV for region 
11, which are consistent with the results using XMM-Newton (the abundance is fixed to 0.1 in 
the fittings of the Chandra spectra). By using the xmmpsf mixing model on the PN spectra 
of the observation on Dec. 8, 2000, we confirm the difference in column density between 
regions 10 and 11. This suggests the hot gas in regions 10 and 11 is likely to be distributed 
in a more or less symmetric manner with respect to the nucleus, and is not confined to the 
galactic disk. Because the near side of the galactic disk is to the SW and the far side to 
the NE, we expect strong absorption by gas in the disk toward region 11 and only weak 
absorption toward region 10, in agreement with our observation. Given the proximity to the 
active jets, it is possible that some of the hot gas in these two regions is heated by the jets. 
This hot gas is also one of the key ingredients in the WYC scenario (see their Figs. 5 and 7), 
in which the gas is ultimately driven into the stellar disk and causes the "line of damage" 
that is seen as the anomalous arms. 

Region 12 includes the northern hot spot (end of the currently active radio jet). The hot 
spot is not resolved from the diffuse emission in the XMM-Newton image. The total emission 
from the hot spot is much lower than that of the more extended emission in this region. The 
point source ~ 10" east of the northern hot spot however, though unresolved, may slightly 
alter the XMM-Newton spectrum from the region. Region 13 locates on a "finger" on the 
SE arm, at a distance from the nucleus similar to region 12. The X-ray "color" suggest 
the temperature in region 12 is likely lower than that in region 13. However, the XMM- 
Newton spectral parameters are almost identical for the two regions. This is likely a result of 
including the hard spectrum point source in region 12. The low absorption column densities 
indicate that the emissions from regions 12 and 13 most likely originate from regions in front 
of the galaxy disk as expected given that the NE side of the galactic disk is the far side. 
Region 14 is includes a "blob" to the east of the nucleus. The Chandra image shows the 
spectrum of the west half of the "blob" is softer than that of the east half. The "blob" 
was reported as a point source (X20) in Vogler & Pietsch (1999). We found the emission is 
resolved with no indication of a point source in the region. The low column density (Table [3D 
to the region suggests the "blob" is located on the near side of the galactic plane. 

Regions 15-18 are chosen to locate on the X-ray "filaments" in the radio "plateau" 
regions. Region 16 also includes the far tip of the northern anomalous arm. The absorption 
and temperature in region 15 is similar to that of region 1. The low absorption indicates the 
emission comes mostly from the near side of the galactic disk. The spectra of regions 16-18, 
on the other hand, show the highest absorption of all regions, ~ 5 — 10 times the Galactic 
value, indicating they are on the far side of the galactic disk. 
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It should be noted that single temperature MEKAL models do not always produce good 
fits to the spectra, as seen in Tabled The fits for the spectra from regions 1, 5, 13, 15, 16, 17 
and 18 are not very good in terms of x 2 P er dof. Most of these regions are in the outskirts of 
the galaxy or in the "plateau" . We introduce a second MEKAL component in these regions, 
assuming that the absorption and abundance of the component are the same as for the first 
MEKAL model. If the second MEKAL is indeed a halo component, the absorption and 
abundance of this component may be different from the first MEKAL component. However, 
the absorption and abundance of the second MEKAL is poorly constrained with the current 
data. Fixing the abundance and absorption does not significantly affect other parameters. 
In Table H] we show the best-fit parameters of the double MEKAL model. By adding a 
component with temperature kT ~ 0.2 — 0.3 keV, the improvement of x 2 is significant. This 
is, of course, the kT ~ 0.22 keV gas found in the integrated spectrum of NGC 4258. The 
normalization of the second MEKAL model is comparable to that of the first component. 
The estimated absorption column densities do not show a significant change from the single 
MEKAL model. 

These results confirm the need for a truly diffuse kT ~ 0.2 — 0.3 keV thermal component. 
This component is observed mostly in the outer part of the galaxy, which is consistent with 
earlier ROSAT results. The nature of this component, however, is unclear . 

As expected, the Chandra spectra of the diffuse regions 12 — 18 agree in general with 
those obtained using XMM-Newton, but are noisier due to the low surface brightness. 

3.3. Hot spots 

We examine the X-ray emission from the radio hot spots. The Gaussian smoothed 
Chandra soft band contours of the northern and southern radio hot spots overlayed on a 
Hubble Space Telescope (HST) WFPC Ha+[NII]A6583+ image are shown in Fig. [6j As 
noticed by Cecil et al. (2000), both radio hot spots are associated with narrow elongated 
enhancements of Ha+[NII] emission, suggestive of shock waves, probably bow shocks driven 
into the ambient gas. The merged Chandra image provides a clear detection of the northern 
hot spot (region 20, Fig. [2]). The centroid of the X-ray emission lies behind (i.e. to the south 
of) the Ha+[NII] emission. The X-ray emission is elongated in the same direction as the 
optical line emission. The hot spot is embedded in a plume of diffuse emission in region 12 
(Fig. [2]). The background subtracted 0.4-2 keV count rate in a 3.5" radius circular region 
centered on the hot spot is (1.1 ± 0.3) x 10 -3 cts/s. The southern hot spot seen in the radio 
and Ha+[NII] images is not detected as a compact X-ray source. However, comparing the 
X-ray and the HST images (Fig. E]), it is easy to see that the X-ray isophotes follow the 
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arc-shaped structure seen in Ha+[NII]. The X-ray emission from the southern bow shock 
also extends to the north to form a shell with a radius of ~ 4.5". The 0.4-2 keV background 
subtracted count rate in this circular region is (2.0 ± 0.4) x 10~ 3 cts/s. We extracted the 
Chandra spectra of the hot spots. While noisy, the spectra show little emission above 2 keV 
in both cases. Because of the small number of counts, spectral modeling is not possible for 
the hot spots. 

A very soft source roughly on the extension of the southern jet (77" or 2.7 kpc south 
of the nucleus, P. A. 176°), previously named as ! SX' in WYC and labeled 19 in Fig. [2J is 
clearly seen in the Chandra and XMM-Newton images, and is speculated to represent another 
jet-ISM interaction. The source is unresolved in these observations. The X-ray spectrum 
using the merged XMM-Newton observations was extracted from a circular region with a 
radius of 20". There is little emission above 2 keV. The spectrum can be well described by 
a bremsstrahlung model absorbed by the Galactic column (allowing the absorption to be 
a free parameter yields a column density consistent with the Galactic value), with kT = 
0.19^08 keV > and /0.2-2 keV = 5.7£g x 10" 15 erg cm" 2 s" 1 . The 0.2-5 keV luminosity would 
be ~ 2.8 x 10 37 ergs -1 if the source is in NGC 4258. 



4. Discussion 

4.1. Geometry of the Anomalous arms 

With spatially resolved XMM-Newton and Chandra spectra, we have shown that the 
absorption to the NW arm is significantly higher than that to the SE arm. This difference in 
absorption has previously been suggested using ROSAT hardness ratios (Cecil et al. 1992). 
A difference in obscuration has also been invoked to explain the observed difference in the 
velocity distribution between the SE and NW arms (Rubin & Graham 1990). Our result 
strongly suggests that the NW arm is partly behind, while the SE arm is in or in front of the 
galactic disk. The simplest interpretation for this is that the arms progressively move out 
of the galaxy disk with increasing distance from the nucleus, as suggested in WYC, where 
the hot gas in the cocoon of the jets hits the galactic disk and shocks the ISM to form the 
observed anomalous arms. In this scenario, the shock heated gas in the outer disk is likely 
to be pushed out of the disk because here the gas is less gravitationally bound than in the 
inner disk. 

The absorbing column to regions 16, 17 and 18 is high, even higher than to the NW 
anomalous arm. One possible explanation for the extra absorption seen to these regions 
is that the gas and dust distributions in the disk are non-uniform, as seen in most spiral 
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galaxies. This is directly testable with IR data. The Spitzer IRAC data has recently become 
available in the archive, which allows us to look at the distribution of warm dust. In Fig. 
we show a color composite of the soft X-ray (Chandra), 8/xm IRAC, and 1.5 GHz high 
resolution radio images. It is clear that the warm dust (traced by the 8/zm IR emission) in 
one of the spiral arms is blocking the soft X-ray emission in the plateau region near the NW 
arm, while the radio emission is not affected. This clearly shows that the northern outer 
arms and the plateau are indeed on the far side of the galactic disk. 



4.2. Hot gas in the anomalous arms 

Our X-ray spectra clearly indicate a thermal nature of the X-ray emission from the 
anomalous arms, and the gas is very likely to be heated by shocks as suggested in earlier 
works. The electron density in the X-ray emitting gas can be estimated using the normaliza- 
tion of the spectra and the volume of the emitting gas. The volumes of spectral extraction 
regions are calculated assuming the regions used for the spectral extractions (Fig. [2]) are 
cylinders with axes along the anomalous arms, except for region 7, for which a spherical 
shape is assumed. With the filling factor / as a free parameter, we estimate the electron 
density N e (filled circles in Fig. [8]) using the PN only normalizations of the single temper- 
ature MEKAL models in Table [31 The density is almost constant (iV e ~ 0.15/" 1 / 2 cm -3 ) 
along the arms. By ignoring the line emissions, we estimate the cooling time of the hot gas, 
which is 

f _ 3N e kT 

i^cool — i.4xio-27 T l/2jvf , 

On the other hand, the evidence of high velocity shock waves found in optical spectra (Cecil 
et al. 2000) allow us to use the best-fit temperature of the hot gas (Table [3]) to estimate 
the velocity of shock waves that heat the gas in the anomalous arms. Assuming the shocks 
are strong (e.g. Hollenbach & McKee 1979), the estimated shock velocity ranges from 500 
to 700 kms -1 (Fig. [8]). The time scale for the shock waves to cross the regions is therefore 
a few times 10 6 yrs, which is much shorter than the cooling time. 

The total volume in our cylindrical regions V ~ 1.2 x 10 65 cm 3 . The total mass of 
the hot gas is therefore ~ 1.5 x 10 7 f 1 ^ 2 M and the total thermal energy is 3N e kTV ~ 
4.5 x 10 55 / 1 / 2 ergs. If the hot gas in the anomalous arms is in thermal equilibrium, the power 
needed to heat the gas should be equal to the radiative loss (assumed to be dominated by 
soft X-ray emission) with L~ 7 x 10 39 ergs" 1 in the anomalous arms. 
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4.3. The nature of the X-ray emission in the hot spots 

To understand the active jets in NGC 4258, it is of critical importance to understand 
the radiation mechanism of the X-ray radiation in the jet hot spots. Unfortunately, without 
good quality X-ray spectra, it is hard to elucidate this mechanism. We briefly discuss the 
possibilities. 

Synchrotron radiation in the X-ray band has been reported in both radio galaxies and 
supernova remnants. In most astrophysical situations, the observed synchrotron spectrum 
is a power-law (S oc v~ a ). Hyman et al. (2001) reported a nonthermal radio spectrum of 
the northern hot spot with flux densities of 0.26 ± 0.08 and 0.71 ± 0.14 mJy at 6 and 20 
cm, respectively. The power-law index was found to be 0.82 ± 0.30. With the Chandra 
count rate, we estimate the radio to X-ray power-law index. Assuming a single power-law 
spectrum with a x = 0.8 (the X-ray flux is only weakly sensitive to a x ), the Chandra count 
rate of the northern hot spot (~ 1.1 x 10~ 3 cts/s) translates to an unabsorbed 0.2-2 keV 
flux of ~ 4 x 10~ 15 erg cm -2 s -1 if the Galactic absorption is assumed, or a flux density of 
2.2 x 10~ 9 Jy at 0.4 keV. Thus the radio to X-ray power-law index is ~ 0.70, which is not 
inconsistent with the radio power-law index, suggesting that synchrotron X-ray radiation 
cannot be ruled out in the hot spots. As noted in § 13. 3[ the X-ray spectra of both radio 
hot spots seem soft, showing little emission above 2 keV. Unfortunately, the current X-ray 
observations do not have sufficient signal-to-noise to exclude an a x ~ 0.7 power-law; longer 
Chandra observations are needed to address this issue. 

The relativistic electrons in the hot spots can up-scatter the cosmic microwave back- 
ground (CMB) photons that dominate the radiation field, into the X-ray band. The relative 
importance of synchrotron and inverse Compton (IC) radiation depends on the relative en- 
ergy densities of the magnetic field and that of the CMB photons. The magnetic field is 
usually estimated with the minimum energy argument, which, to the order of unity, gives 
equipartition of energy between magnetic field and the relativistic particles. If we assume 
cut-off frequencies of 10 7 and 10 11 Hz, adopt the spectral parameters from Hyman et al 
(2001), and choose a radius of the radiation region of the northern hot spot of 2" based 
on the high resolution radio image from Cecil et al. (2000), we find an equipartition field 
B m w 21 /iG (Ginsburg & Syrovatski 1964), which is ~ 7 times the equivalent magnetic field 
of the CMB. Therefore, it is unlikely that the X-ray radiation is dominated by IC, unless the 
magnetic field in the hot spots is at least an order of magnitude lower than B rn . The same 
argument applies to the southern hot spot. 

Optical spectra of the emission associated with the radio hot spots (Cecil et al. 2000) 
suggest the gas is heated by shocks with velocity of 350 ± 100 kins -1 . Such shocks are 
capable of heating gas to X-ray emitting temperatures. This seems to agree with the X- 
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ray spectra if the temperature of the hot gas is below 2 keV. The clear spatial correlation 
between the X-ray and optical emission (Fig. [6]) suggests that the X-ray emission is likely to 
be associated with these shocks. If the X-ray emission comes from shock heated plasma, and 
could be modeled as bremsstrahlung absorbed by the Galactic column with a temperature 
of 0.5 keV (the X-ray colors of the hot spots appear similar to that of the anomalous arms), 
the count rates would translate to fluxes of ~ 3.9 x 10 -15 and ~ 7.0 x 10 -15 erg cm -2 s -1 , or 
absorption corrected luminosities of 2.8 x 10 37 and 4.7 x 10 37 ergs -1 between 0.4 and 2 keV, 
for the northern and southern hot spots respectively. These X-ray luminosities are more 
than a factor of 100 below that of the anomalous arms (§4.2). 



4.4. Jet energetics and the anomalous arms 

If the X-ray emission in the hot spots is confirmed to be thermal in nature, the emission is 
likely to come from the post-shock gas heated by the bow shocks that propagate into the ISM. 
We favor this interpretation because: 1) The X-ray emission show excellent agreement with 
the shock structures identified in HST images; 2) The X-ray color suggests the temperature 
of the shock heated gas is ~ 0.5 keV, which implies a shock velocity of 650 km s -1 (e.g. 
Hollenbach & Mackee 1979) and which is compatible with the optically measured range of 
shock velocities (350 ± 100 kms -1 , Cecil et al. 2000). The total power in the shocks can be 
estimated as follows. Using the X-ray fluxes in the hot spots, we found the electron density 
in the post-shock gas to be 0.20/ -1//2 and 0.19/ -1 / 2 cm -2 ( / is the filling factor) for the N 
and S hot spots, respectively (we have adopted spherical volumes with radii of 3.5"=122 pc 
and 4.5"=157 pc for the northern and southern hot spots, respectively). The total power of 
the bow shocks would be |/w 3 v4, where p is the density of the ISM, which is 1/4 the density 
of the shock heated gas in the hot spots in a strong shock, v s is the velocity of the shock, 
and A is the area of the hot spot. Assuming v s = 350 kms -1 , the power of the bow shocks 
is 8.5 x 10 38 /~ 1//2 and 1.7 x 10 39 / -1//2 ergs -1 for N and S hot spots, respectively. The total 
power is therefore ~ 2.6 x 10 39 / -1 / 2 ergs -1 . This power is much higher than the X-ray 
luminosity of the hot spots (7.5 x 10 37 ergs -1 , see § 14.31) . In the case of light, non-relativistic 
jets, the shock heated ISM mixes with the back flow from the post shock jet gas (Williams 
1991). The gas then fills a cocoon around the jets which may stabilize the jet flow. Almost all 
the power of the jets is transported to the cocoon. It is interesting to note that the power in 
the bow shocks driven by the jets (~ 2.6 x 10 39 / -1 / 2 erg s -1 is of the same order as magnitude 
of the power required to heat the hot gas in the anomalous arms (~ 7 x 10 39 ergs -1 ; § 14. 2p . 
This is encouraging for the WYC model, in which the anomalous arms are powered by the 
shock waves driven into the galactic disk by the out-of-plane jets. The shock waves in this 
scenario could be a result of the cocoon expanding into the ISM; the shock waves transport 
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part of the jet energy in directions transverse to the jet direction. In particular, the shock 
waves caused by the expanding cocoon into the disk gas may be responsible for the arc-like 
features seen in the Chandra image (regions 4 and 7 in Fig. [2]). These arc-like features are 
at a similar nuclear distance to the jet hot spots N and SX (regions 20 and 19) and may 
represent the impact of "splatter" from the ends of the jets onto the galactic disk. 

4.5. Jets and the accretion disk 

If the currently active jets are indeed responsible for heating the anomalous arms, in 
particular the inner arms in the WYC scenario, we argue that the power of the jets should 
be comparable to or greater than the X-ray luminosity of the anomalous arms. This is 
particularly interesting in NGC 4258 because the accretion power in the AGN can be directly 
inferred from the disk structure inferred from water maser observations. This allows us to 
constrain the fractional accretion power which has been injected into the radio jets, which 
has important implications for the physics of accretion disk. 

Assuming a thin warped disk illuminated by a central X-ray source (Miyoshi et al. 
1995; Greenhill et al. 1995), and using the masing conditions, Neufeld & Maloney (1995) 
found an accretion rate M = 7 x 10~ 5 a M yr -1 , where a < 1 is the usual dimensionless 
viscosity defined in Shakura & Sunyaev (1973). The accretion power is therefore P = 
r]Mc 2 = 4 x 10 41 a erg s , where we have assumed a standard thin accretion disk with 
efficiency 77 = 0.1. On the other hand, the total X-ray luminosity of the inner anomalous 
arms (regions 4, 5, 6, 7) is ~ 4 x 10 39 ergs -1 . We further assume that only ~ 50% of 
the energy from the jets has been converted to thermal energy in the disk because only jet 
power driven down towards the disk can heat the anomalous arms. The power of the jets is 
thus > 8 x 10 39 ergs -1 . This result suggests that the power in the jets is > 2a _1 % of the 
accretion power. The power of the jets thus estimated is very similar to the X-ray luminosity 
of the nucleus. Note that the argument above also applies to models other than the WYC 
scenario as long as the jets from the AGN are the sources of the thermal energy in anomalous 
arms. The only difference is the factor of 2 introduced above to account for the out-of-plane 
jets. 

Using a sample of 9 elliptical galaxies, Allen et al. (2006) have found recently a tight 
correlation between accretion and power of the jets (inferred from X-ray emitting gas in the 
halos/bubbles in the galaxies), 

log(P/10 43 erg s -1 ) = A + B log(P jct /10 43 erg s -1 ) (2) 

where P = t]Mc 2 , rj = 0.1, A = 0.62±0.15 and B = 0.77±0.18. Since the powers of the jets 
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in this sample are typically in the neighborhood of 10 43 ergs -1 , it is interesting to test if 
the relation extends to low luminosity AGNs. In NGC 4258, we found log(P/10 43 ergs -1 ) = 
—1.4 + log a. On the other hand, using the power of the jet estimated above, we obtain the 
RHS > —1.76. It is obvious that Eqn. [2] is satisfied in NGC 4258 to within the uncertainty 
of the relation. Should our estimates of the jet power be confirmed in future observations, 
this result would suggest: 1) While the accretion mode may be very different in elliptical 
galaxies and low luminosity AGNs, like NGC 4258, the physics of jet formation is likely to 
be the same; 2) Since the viscous time scale at 0.13 pc is long (~ 10 8 yr if a = 0.1, Gammie, 
Narayan, & Blandford 1999), the accretion rate must have been stable over such a long time 
scale. 



5. Concluding remarks 

In this paper we have examined the spatially resolved archival Chandra and XMM- 
Newton spectra of NGC 4258 in detail. We summarize our main results. 

1. Using stacked XMM-Newton spectra, we present a time averaged spectrum of the 
nucleus of NGC 4258 over 1.5 yrs. Our spectral parameters agree very well with those 
reported in Hanson06, who use a slightly different method. We confirm the detection 
of a weak Fe Ka line in the combined spectrum, with En ne = 6.39 ± 0.04 keV, and 
EW = 40 ± 33 eV. 

2. The Chandra spectra of the off-nuclear source J121857.3+471812 show that the source 
is responsible for the second power-law component detected in the XMM-Newton spec- 
trum of the nucleus. 

3. The X-ray spectra of the anomalous arms confirm the thermal nature of the X-ray 
emission. The temperature of the hot gas (kT ~ 0.5 keV) agrees with the idea that 
the gas is heated by shocks with velocity similar to those implied by optical kinematic 
measurements (Cecil et al. 1995). 

4. The absorption column density to the SE anomalous arm is consistent with the Galactic 
value, while the column density to the NW arm is significantly higher than the Galactic 
value. This suggests the SE arm is in the galaxy disk or slightly in front of it, and, 
part of the NW arm is on the far side of the galactic disk of NGC 4258. The X-ray 
emission associated with the NW plateau region appears to be behind the galaxy disk, 
as inferred from the X-ray column densities and Spitzer IRAC 8/xm images. 
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5. The absorption column density to the diffuse region SW of the nucleus is higher than 
that to the diffuse region NE of the nucleus, as expected from the fact that the SW 
side of the galactic disk is the near side. This indicates that a significant part of the 
hot gas within ~ 2 kpc of the nucleus is out of the disk plane. 

6. We confirm that a kT ~ 0.2 — 0.3 keV diffuse thermal component is present in some 
regions of the outer part of the galaxy. 

7. We estimate the number density of the X-ray emitting gas along the anomalous arms 
is ~ 0.15/ -1 / 2 cm" 3 (/ is the filling factor) and the velocity of shock waves that heats 
the anomalous arms is 600 — 700 kms" 1 . The cooling time of the hot gas is much 
longer than the dynamical time scale of the shock waves. 

8. The Chandra image of the radio hot spots shows that the X-ray emission is closely 
associated with the shock structures identified in HST images. While synchrotron 
radiation cannot be ruled out, it is most likely that the X-ray emission in the hot spots 
is thermal in nature, and is heated by the shocks responsible for the emission lines seen 
at optical wavelengths. 

9. Assuming the X-ray emissions from the radio hot spots at the ends of the radio jets 
come from shock heated gas, we found the power in the shocks is comparable to the 
power needed to heat the anomalous arms, suggesting that the currently active jets 
may be responsible for heating the gas in the anomalous arms. If the AGN jets are 
indeed the source of thermal energy in the anomalous arms, we estimate the power 
of the jets to be > 8 x 10 39 ergs -1 . The power of the jets estimated in this way is 
> 2a _1 % (where a is the viscosity parameter for the disk) of the accretion power in 
the disk found by Neufeld & Maloney (1995). This seems to agree with an accretion 
power versus jet power relation found in elliptical galaxies, even though NGC 4258 has 
an AGN luminosity three orders of magnitude lower than studied in this relationship. 

We have shown that the currently active jet could be responsible for the formation of 
the anomalous arms. The X-ray emission from the anomalous arms therefore provides a 
useful perspective on the energetics of the radio jets and their relation to the accretion disk 
in NGC 4258. Our results are consistent with the scenario proposed in WYC, which suggests 
the anomalous arms are caused by the impact of shock waves from the jet driven gas in the 
halo onto the galactic disk. The arc-like structures in the anomalous arms (Fig. [Tj) seen in 
the Chandra image could be a natural consequence of hot halo gas shocking the ISM in the 
galaxy disk. 

To make further progress, good quality X-ray spectra of the radio hot spots are needed. 
The spectra will be useful to determine the radiation mechanism, and hence gain a better 
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understanding of the jets. Numerical simulations of highly inclined radio jets in and above 
a dense galaxy disk would be desirable to better quantify the energy transport. 

We thank Dr. Keith Arnaud for advice on XSPEC. This work is supported by NASA 
through LTSA grant NAG 513065 to the University of Maryland. 
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Table 1. Observations of NGC 4258 



Observatory Obs ID 


Detector (filter) a 


Observation Starting Time 


Exposure time (ks) b 


Chandra 350 


ACIS-S3 


2000-04-17 10:52:45 


14.04 


1618 


ACIS-S3 


2001-05-28 21:21:44 


20.93 


XMM-Newton 0059140101 


MOS1 (Medium) 


2001-05-06 18:56:45 


12.00 




MOS2 (Medium) 


2001-05-06 18:56:42 


12.00 




PN (Medium) 


2001-05-06 19:36:07 


7.87 


0059140201 


MOS1 (Medium) 


2001-06-17 16:05:37 


12.44 




MOS2 (Medium) 


2001-06-17 16:05:37 


12.44 




PN (Medium) 


2001-06-17 16:43:12 


9.05 


0059140401 


MOS1 (Medium) 


2001-12-17 05:26:41 


6.32 




MOS2 (Medium) 


2001-12-17 05:26:41 


6.83 




PN (Medium) 


2001-12-17 06:05:54 


3.88 


0059140901 


MOS1 (Medium) 


2002-05-22 21:35:06 


14.86 




MOS2 (Medium) 


2002-05-22 21:35:06 


14.96 




PN (Medium) 


2002-05-22 22:08:22 


10.06 


0110920101 


MOS1 (Medium) 


2000-12-08 22:25:31 


20.24 




MOS2 (Medium) 


2000-12-08 22:25:31 


20.26 




PN (Third) 


2000-12-08 23:32:01 


14.44 



a The CCD on which NGC 4258 is observed. 
b High background intervals have been filtered. 
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Table 2. Model parameters for the X-ray spectra of the nuclear region from it 

XMM-Newton observations 



Darameter name 


FIT 1 FIT 2 


N Htgal (10 20 cm- 2 ) 


2 2Q 2 - 70 2 36 3 - 36 


kT (keV) 
Abundance 

K a 


o.578;Il o.57g:ig 
o.i3g:ig o.io^l 

3.89|:|1 x 10" 4 4.08f;|f x 10" 4 


N m (10 22 cm" 2 ) 
Ti 


i n 1 10.3 o 7o9.26 
1U.1 848 8 -' 8 8.31 

1.96 2 'i 1.66 4 5 1 

2.80|i| x io- 3 1.82 2 ;!! x 10- 3 


jV# 2 (10 22 cm" 2 ) 

r 2 


0.34 (frozen) 
0.10^ 4 20 1.74 (frozen) 
1.34 4 :^ x 10~ 5 3.36||| x 1(T 5 


E Fe (keV) 
a (keV) 

fline 


6 3 7 6 - 42 fi 37 6 - 42 

u -°'6.30 u -°'6.31 

n nc;0.i5 q 040.14 
u - uo o.oo u - u ^o.oo 

3-50f;|| x 1(T 6 1.82?;{g x 1(T 6 


X 2 /dof 


856/836 863/837 



& K = 10- 14 /(47r( J D A (l + z)f) J n e n H dV, where D A 
is the angular size distance to the source (cm), n e is the 
electron density (cm -3 ), and n# is the hydrogen density 

(cm" 3 ) 

b defined as the count rate at 1 keV 
c unit: photons cm -2 s _1 



Table 3. Model parameters for the X-ray spectra of extended regions: Single temperature MEKAL models 



PN PN + MOS 



Region No. 




kT (kcV) 


abundance 


K (xl0 -4 ) b 


X 2 /dof 




kT (keV) 


abundance 


K (xl0^ 4 ) b 


X 2 /dof 


1 


a oil 68 
U -81o:28 


o.4ig:$ 


0.068:8! 


3-44 4 i 7 


97.3/91 


0.008:88 


a o 4 ^3 

0.53 : 47 


„ „n f!7 

0-068:81 


z - ' °2.62 


200/163 


2 


1.46g;f§ 


0-428;|l 


0.058:81 


2-3o 2 :8t 


56.1/68 


o.4oJ:88 


o.498:i 


0-068:81 


l-85 2 i8 


114/122 


3 


1 1 1 l - ' 4 

J -- ± - L 0.60 


0.49g : l 7 




4.59S:& 


157/176 


0-31o:oo 


a rc*i 52 

0-51o!50 


o.n8:io 


3-76 4 i8 


287/258 


4 




0.49g : l 7 


o-i28:io 


°. 5,3 3.36 


162/180 


a ooO 77 

"••J^oioo 


A c 1 52 


o.i48:l| 


1 o3 45 
• : >.l z 2:88 


275/261 


5 


i ci 2.05 
1 - oi 1.02 


A coO.bd 
U - OZ 0.50 


o.i38:i8 


5-88i:i 


248/209 


o.99j : |8 


n cq0.55 
u - oo 0.52 


o.i28:l? 


5.32 4 : 97 


374/301 


6 C 


3.05 2 ;7q 


o-568;|l 


o-io8:o9 


i4-2li:I 


992/986 


O QC3.20 

Z.«D 2 . 51 


0.5781 


0108:^9 


13.7118 


1591/1481 


7 


4.696;38 


o-548:|8 




i.74?:£ 


51.8/86 


^• 00 2.96 


o.538:i8 


0.118:09 


i-66?:§| 


88.7/139 


8 


3 98 4 ?g 






5 42§ ?§ 

J -^ z 4.17 


1 94/1 96 


3 313-95 


52S-5 3 . 

u - oz 0.50 


1324^ 

U . 1J 0.12 


4 725 i? 

' z 4.25 


332/286 


9 




o.498:I 


o.o98: I 


3-24|:i 


100/113 


4-36|:|° 


a ci 0.54 
u . oi 0.49 


o.i28:i5 


2.48 2 : ?f 


205/181 


10 c 


2.45?:|g 


n MO-58 
u - 04 0.54 


o-098: g 


i2.8l?:l 


723/853 


2.04f-3| 


o.578:i 


o.o98: § 


1.18I:?! 


1323/1339 


ll c 


5.24|it 


n cc0.58 
u - oo 0.55 


o.o98: 8 


q 1 5 104 

y - 1,J 8.22 


553/571 


4.9it:I^ 


0-588:1? 


10° 11 

U - 1U 0.09 


078.71 

O.d / 7 4 8 


903/928 


12 


1 69 2 52 


u - 01 0.51 


o.iog:Si 


q 7q4.29 
°- ' °3.26 


126/156 


J^- ' 1.08 


a c ^0.56 
u -°%.52 


10° 11 

U - 1U 0.09 


3 .6 27 


212/238 


13 


i cn2.03 
1.0\) 1 19 


fl fSI 0.54 
U - Oi 0.51 


0.098:0s 1 


8-54?:i 


261/224 


i.05 : 3 7 I 


a cq0.54 
U - O,3 0.52 


n in 11 

U - 1U 0.09 


7.66?i2 


461/325 


14 


u - oo 0.07 


o.528:iB 


1 1 013 
u - ii o.09 


9 c;q3.i6 

z - oy 2.29 


124/149 


0.84oJ 4 


a C4O.59 
u -°%.53 


o.ii8:18 


9 CQ2.87 
z .°°2.19 


225/226 


15 


n 47 101 

0.00 


o-468i| 


0-078:81 


4.64l;l| 


202/168 


n nn - 36 
u - uu o.oo 


o.5o8iI 


0.088:8? 


A 1Q4.53 

4.iy 4 04 


339/251 


16 


o 91 10.0 
°- zi 7.20 


o.488:tI 


0.088:88 




284/208 


6-211:1 


o.598:f 7 


0.108:^ 


7 n7 7 - 78 

'■ u '6.42 


466/312 


17 


6.40 5 7 ;?| 


0-508:1^ 


o.io8: i 


6.i6l:? 2 7 


211/193 


4.49l: 2 7 l 


o.588:f 7 


0.130^2 


4.57|:I 4 


345/286 


18 


10.9 7 4 7 2 


o.48g:i 


0-068:81 


7.45181 


183/164 


7.74« : ? 2 


o.6o8:i 


o.iogii 


4.63i:?| 


326/259 



a in units 10 20 cm 2 

h K = 10 _14 /(47r(I? y i(l + z)) 2 ) J n e riHdV, where Da is the angular size distance to the source (cm), n e is the electron density (cm -3 ), 
and nn is the hydrogen density (cm -3 ) 

c spectra in 0.2-10 keV band were fitted (see text for detail). 
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Table 4. Model parameters for the X-ray spectra of extended regions: Two temperature 

MEKAL models 



Region No. N H a kT (keV) abundance K (x!CT 4 ) b kT (keV) K (x!0~ 4 ) b x 2 /do/ 



PN 



1 


S5 1 - 60 

u - od 0.00 


U - OD 0.49 


0. 


1 90.20 

lz 0.10 


1 4.D 194 
1 - 4U 1.03 


n iq - 23 

u - iy 0.12 


-1 9^1.65 
1 - zd 0.94 


77.4/89 


5 


i oqi.88 


0.56°;i 


0. 


170.21 
1 '0.14 


3.87 4 -i3 


fl 910.26 
u - zl 0.16 


1 141-65 


221/206 


13 


1 - z ' 1 0.87 


0-59°;!? 


0. 


1K0.17 
lo 0.13 


4.89i;lf 


9^0.26 
u - zo 0.20 


9 Ifi 2 - 72 
z - iD 1.76 


202/222 


15 


0.43^ 


U.OSq 54 


0. 


iqo.27 

ly 0.15 


i-87?:i| 


170-20 
u - i '0.14 


1 70 2 - 27 

L - ' u 1.31 


132/166 


16 


9-48^1 


U - D4 0.60 


0. 


1 7O.22 
i '0.14 


4.42|;i 


0.17»; 2 i 


4 9n 609 

^• zu 3.06 


188/206 


17 


c ni7.61 
D.U1 4 64 


fl fil 064 
U - O1 0.58 


0. 


1 R0.23 
iS 0.14 


3.24 4 ;°| 


0.18g; 2 i 




174/191 


18 


11 4 14 - 2 
1J - 4 8.60 


770-86 
u - ' '0.62 


0. 


1 0O.I8 
lo 0.09 


2-50f; 4 | 


0.30g;| 4 


4 176.20 
' 2.66 


144/162 



PN+MOS 



1 


n nn - 85 

u - uu 0.00 


n ^70-60 

u -°'0.55 


0. 


190.16 

lz 0.10 


1 -«l-47 


n 19 - 22 
u.iy .i4 


1 19 131 
1 - lz 0.98 


163/161 


5 


Q2 1 - 29 

u - yz 0.53 


0-57°;i 


0. 


1 00.20 

i8 0.16 


3.70^ 24 


20°- 24 

u - zu 0.18 


i-025:ii 


331/299 


13 


0-92 ;| 7 


0-590:1° 


0. 


1R0.17 
iD 0.14 


4.9lfct? 


n 99024 

u - zz 0.20 


i-82 2 :M 


366/323 


15 


0-00° : o 7 


0-59g;i 


0. 


iqo.22 

ly 0.16 


1.86 2 i| 


0-188;?g 


L4 1 : l 


220/249 


16 


9-16 7 ° 97 


0-63g : gS 


0. 


iD 0.14 


4-82|; 4 | 


0-228$ 


3.452 77 


334/310 


17 


5 SI 6.67 


R10.64 
U - O1 0.59 


0. 


1 sO.22 
iS 0.16 


Q13.47 
°- oi 2.68 


94°-26 

u - z ^0.20 


1.41^0 


279/284 


18 


iu - u 8.29 


n 770-82 
u - ' '0.73 


0. 


1 cO.19 
io 0.12 


9 91 2.63 
z - zi 1.80 


A 090.35 

u - oz 0.30 


3-06|: 4 | 


261/257 



a in units 10 20 cm 2 

h K = 10 _14 /(47t(.Da(1 + z )) 2 ) J n e riHdV, where Da is the angular size distance to the source 
(cm), n e is the electron density (cm -3 ), and ur is the hydrogen density (cto~ 3 ) 
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Fig. 1. — The adaptively smoothed Chandra true color image of NGC 4258. The RGB colors 
are assigned to represent soft X-ray bands defined as: red: 0.4-0.7 keV; Green: 0.7-1.4 keV; 
Blue: 1.4-2 keV. The northern hot spot at the end of the radio jet is marked N. 
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Fig. 2.— Adaptively smoothed Chandra 0.4-2 keV (left) and XMM-Newton 0.2-2 keV (right) 
images. The regions for spectral extraction are also displayed. 
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PN-Black; MOS-Red 




channel energy (keV) 

Fig. 3. — The stacked 2.2-10 keV PN and MOS spectra and the best-fit model of the nuclear 
region. The spectra are rebinned so that the signal-to-noise ratio in each bin is > 9 in 
XSPEC for better visualization. An Fe Ka emission line is clearly visible in both spectra. 
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Off nuclear source 



o 
o 




Energy (keV) 

Fig. 4. — The Chandra spectra of the off-nuclear source. Red: spectrum using the May 28, 
2001 observation; Black: spectrum using the April 17, 2000 observation. The spectra are 
grouped so that each bin has > 3 cts. A 6.4 keV Fe Ka emission line is added to the model. 
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Fig. 5. — The absorption column densities and temperatures along the anomalous arms. 
Region numbers are defined in Fig. [2J The empty and filled circles represent respectively 
the best-fit using PN spectra only, and both PN and MOS spectra. The small diamonds 
represent the best-fits using combined Chandra spectra. The column density is higher to the 
NW anomalous arm than to the SE anomalous arm, indicating that the former is on the far 
side of the galactic disk. 



Fig. 6. — The Chandra contours of the radio hot spots (left: northern hot spot, marked N in 
Fig. 1; right: southern hot spot) overlayed on the HST Ha+[NII] image (Cecil et al. 2000). 
The southern hot spot is marked S in Fig. 5 of WYC. A clear correspondence between X-ray 
and optical emissions is seen in these putative shocks. The high peaks in the X-ray contours 
that correspond to the optically identified shocks are marked "H" , while the "valley" behind 
(i.e. to the north of) the S shock and other areas of low X-ray emission around the S shock 
are marked with "L" . 
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Fig. 7. — The composite VLA 1.46 GHz (Cecil et al. 2000) (blue), Chandra (red) and Spitzer 
8//m (green) image of NGC 4258. The X-ray emission at the northern end of the anomalous 
arms and northern plateau regions is clearly obscured by the warm dust in the galactic disk 
shown in the 8/zm image. 
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500 F. H 

SE1 23456789 NW 

Region Number 

Fig. 8. — The density of the hot gas in the anomalous arms (upper panel) and the derived 
shock velocity using the measured X-ray temperatures, assuming the gas is heated by strong 
shocks (lower panel) . The dashed line is the mean value of NJ 1 ' 2 . 



